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We study the finite-temperature phase of a gluon ensemble in a variational approximation to QCD 
in the Coulomb gauge. We derive and numerically solve the underlying Dyson-Schwinger equations 
up to one-loop order. Assuming the subcritical solution at T = 0, we find a sharp transition in the 
infrared value of the gluon energy at a critical temperature. 
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I. INTRODUCTION 

The determination of the phases of hadronic matter 
plays a major role in understanding the mechanisms of 
confinement and dynamical symmetry breaking in Quan- 
tum Chromodynamics (QCD). The methods available 
to investigate these phase transitions are the following: 
There are lattice simulations, which have been successful 
in mapping out the deconfinement transition of the QCD 
phase diagram as a function of temperature for near zero 
chemical potential [1-5], and phenomenological models, 
that in addition can cover the high-density regime [6- 
10]. Finally, in the asymptotically large temperature or 
density limit, due to asymptotic freedom, the weak in- 
teractions between quarks and gluons are expected to 
determine the properties of the quark-gluon plasma [11- 

In this paper we investigate the thermal properties of 
the low-density phase using a set of tools that bridge 
QCD and phenomenology. In particular, we formulate 
the problem in the physical. Coulomb gauge, canoni- 
cal Hamiltonian framework of the pure gauge theory. 
While there have been numerous studies of QCD based on 
Dyson-Schwinger resummation techniques [16-19], renor- 
malization group flow equations [20] , and lattice simula- 
tions [21-23] in covariant gauges, the few that exist in 
physical gauges are rather loosely related to the underly- 
ing QCD interactions [24-29], with a recent attempt at 
a self-consistent calculation at finite density [30].^ 

The advantages of physical gauges for phenomenology 
and for developing physical intuition are clear, and we 
summarize them here. The degrees of freedom of the pure 
Yang-Mills (YM) theory are transverse gluons, and ther- 
mal excitations connect color singlet states of arbitrary 
number of gluons. Transverse gluons are expected to be 
effective only at high temperatures, while at low temper- 
atures it would be more effective to compute the partition 
function in terms of the ground state glueballs [34, 35] . 
The underlying interactions in Coulomb gauge are dom- 
inated by the instantaneous Coulomb potential acting 
between color charges. In the non-Abelian theory, the 



potential not only couples charges but it also depends 
on the gluon distribution of the state in which it is cal- 
culated. At zero temperature, in the vacuum state this 
distribution is such that the Coulomb potential becomes 
confining, i.e., proportional to the distance R between 
the external color charges, V{R) = acR [36, 37]. Using 
various approximate, variational models for the ground 
state YM wave functional, it has been possible to ob- 
tain a potential which is confining [38] or almost con- 
fining, i.e., V{R) ex R^-' with e w 0(10%) [39-41]. 
The Coulomb string tension ac is larger than the string 
tension computed from the temporal Wilson loop. This 
is because the Coulomb potential represents the energy 
of a static quark-antiquark pair submersed in the QCD 
vacuum, while the Wilson loop measures the energy of 
the exact QQ state in which the gluon distribution is 
squeezed by closed vortex lines. Since the Coulomb po- 
tential is an instantaneous observable, one might expect 
that it remains confining even in the high-temperature 
limit [37]: At high temperatures the integration over 
transverse fields becomes even less restricted than in the 
vacuum, and, according to the Gribov-Zwanziger confine- 
ment scenario [42, 43], Coulomb confinement originates 
from large field configurations near the the Gribov hori- 
zon. 



^ At zero temperature and density, Dyson-Schwinger studies in 
Coulomb gauge have been performed in Refs. [31-33]. 



In the following, we investigate the finite temperature 
properties of Coulomb gauge Yang-Mills theory with fo- 
cus on the aspects of deconfinement at finite temperature. 
We extend the variational approach of Refs. [38, 40, 41] to 
finite temperature. In particular, the variational Gaus- 
sian ansatz for the vacuum wave functional is extended 
to include single particle, quasi-gluon excitations. In Sec- 
tions II and III we present the general setting for the fi- 
nite temperature, canonical Coulomb gauge problem. In 
Section IV we discuss the details of the variational ap- 
proximation. In Section V we give details of the numeri- 
cal computations and results. Our summary and outlook 
are given in Section VI. 
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II. HAMILTONIAN APPROACH AT FINITE 
TEMPERATURES 

After resolving Gauss's law in Coulomb gauge, the 
Yang-Mills Hamiltonian reads 

Hru - ^ /d'a; {J-'[A] nj[A] n + B2) + He 



Hk + He 



He. 



He = ^Jd^xd^y J-'[A] p"(x) J[A] <(x,y) p^(y), 

(2) 
where n"(x) = —iS/SA"'{x) is the canonical momentum 
(electric field) operator, and 

J[A] = Det(-DV) (3) 

is the Faddeev-Popov determinant with 

D = V + .9A, A"'' = fcA", (fc)"^ = r"'' (4) 

being the covariant derivative in the adjoint representa- 
tion. Furthermore, 



p''{x)^-r'"'A'' -w 



(5) 



is the color charge density of the gluons and 

FlV,y) = (x,a|(-DV)-i(-V2)(-DV)-i|y,&> (6) 

is the so-called Coulomb kernel. Its vacuum expectation 
value (i^^''(x, y)) represents the static non-Abelian color 
Coulomb potential. 

The gauge fixed Hamiltonian Eq. (1) is highly non- 
local due to Coulomb kernel Fa(x, y), Eq. (6), and due 
to the Faddeev-Popov determinant J[A], Eq. (3). In ad- 
dition, the latter also occurs in the functional integration 
measure of the scalar product of the Coulomb gauge wave 
functionals 



{i^l\OH2 



DAJ[A]iPl[A]0'ip2[A]. 



(7) 



In Ref. [40] the Yang-Mills Schrodinger equation was 
solved by the variational principle using the following 
ansatz for the vacuum wave functional 



(A|0) 



1 



VJ[a] 



im, 



{A\0) ==A/'exp ( -^ / dk A{-k)uj{k)A(k) ] , 



(8) 



where 



dk 



d^k 

(2^ 



(9) 



The pre-exponential factor removes the Faddev-Popov 
determinant from the scalar product Eq. (7). The kernel 



Lu(\i.) was determined by minimizing the vacuum energy 
(Hym). which yields an a;(k) which can be well fitted by 
Gribov's formula 



c.(k) 



M* 



(10) 



and which is in satisfactory agreement with the lattice 
data [44], with M « 860 MeV. 

The present paper is devoted to study Yang-Mills the- 
ory at finite temperatures, which is defined by the density 
operator 

V^Z-^cxpi-l3Hy,,), (11) 

where /3 ~ \/T is the inverse temperature and 

Z^Tre-'^^™ (12) 

is the partition function. 

To calculate the trace in the thermal averages 



(O) = Tr(02?) 



(13) 



we need a suitable basis in the gluonic Fock space, which 
we choose as follows: Let af(k) be the operator which 
annihilates the vacuum state |0) [Eq. (8)], i.e.. 



<(k)|6) = o. 



(14) 



Then a complete basis in the gluonic Fock space is given 

by 



{|n)} = {|0),af(k)|0), ar(k)af(q)|0), ...}. (15) 

Following Ref. [40] we choose the basis states of the glu- 
onic Fock space in the form (cf. Eq. (8)) 



{|n)} = {j-V2[A]|n)}. 



(16) 



The thermal expectation value Eq. (13) can then be ex- 
pressed as 



(0)=tr(Pd), 
where the operation '^' is defined by 

= ji/2[A]OJ-i/2[A] 



(17) 



(18) 



and 'Tr' means that the trace is evaluated in the basis 
of the states {[«)} [Eq. (15)]. The transformed density 
operator reads explicitly 



V = Z-^ cxp{-/3Hym), Z = tre^'^-f^™. 



(19) 



This operator is too difficult to handle in semi-analytical 
calculations. In analogy to the zero-temperature case, 
where the Gaussian vacuum wave functional Eq. (8) was 
assumed, we will replace the exact (transformed accord- 
ing to Eq. (18)) Yang-Mills Hamiltonian Hym by a single- 
particle operator 



bt/ 



h^ dkn{k)a'P{k)atik), 



(20) 



where the kernel J7(k) will be determined by minimizing 
the free energy 

T = (i^YM) - TS. (21) 

Here S is the entropy, which is defined by 

5'^-trPlnP. (22) 

By straightforward manipulations the expression for the 
entropy Eq. (22) can be cast into the form 



S' = lnZ-;3 



ainZ 



(23) 



which will be convenient in later calculations. 

III. COLOR PROJECTION 

A. Exact projection 

By definition the trace in the thermal averages Eq. (17) 
should be taken in the physical Hilbert space. Before 
gauge fixing the physical Hilbert space is given by all 
gauge invariant states. After resolving Gauss' law in 
Coulomb gauge, the physical Hilbert space is defined by a 
complete set of wave functionals of the transversal gauge 
field that are invariant under global gauge transforma- 
tions (the latter is not fixed by the Coulomb gauge con- 
dition). These states are annihilated by the total color 
charge operator^. 



Q" 



A^xp^i-K). 



(24) 



However, an individual basis state of the set Eq. (15) 
will, in general, carry a non-zero color charge, and the 
use of the basis Eq. (16) will lead to a colored statistical 
ensemble. Therefore we project these states onto color 
singlet states using the projector 



V= I Ayi{e) exp {idaQ''\ , 



(25) 



where dn{6) denotes the Haar measure of the gauge 
group parametrized in terms of the color angles 0a ■ The 
thermal average projected onto zero-color states reads 



(O) =fr(OPP), 



(26) 



where V is given by Eq. (19) with Hyu replaced by h 

[Eq. (20)]: 



- y-^c-Ph 



V = Z 



Z = Tt{c-^'^V) 



^ Also in the functional integral formulation after fixing to 
Coulomb gauge a careful treatment of the zero modes of the 
Faddeev-Popov operator related to the global gauge transforma- 
tions constrains the ensemble of transversal gauge fields to those 
with vanishing total color charge [45]. 



The density operator V [Eq. (27)] is color singlet and 
hence commutes with the total color charge operator Q" 
[Eq. (24)], which in terms of the creation and annihilation 
operators reads 



Q-^ir"^ dfca^T (k)<(k) 



bt^ 



(28) 



With the explicit form of the projector Eq. (25) we have 



VV^ f dfi{e)Ve, 



where 



Va 



js- 



Qf) 



-Ph+iO-Q 



(29) 



(30) 



is the density operator in the presence of an external color 
field (— i6'aQ"//3), i.e., for fixed color angle 6a- Due to 
the presence of the external color field (— i6'aQ"//3), this 
density matrix is non-diagonal in color space. However, 
since the total charge operator Qa is hermitian and h is 
color singlet, we can diagonalize Ve. For simplicity, we 
consider the gauge group SU{2). Then wc may write 



Vg^uHe)VeU{e), 



(31) 



where U{e) lives in the coset SU{2)/U{1) and 

Vg ^cxp{- f3h + i9Q^) (32) 

lives in the Abelian subgroup. In a parametrization of 
the gauge group SU{2) corresponding to the coset de- 
composition Eq. (31) the Haar measure reads 



dp{e)= fdfi(e) f dp{e), 
MO) = ■ 



Is- 
d9 sin^ 



(33) 



2' 



where d/z(0) denotes the measure for the integration over 
the coset's SU{2)/U(1) ~ S^ degrees of freedom. 

In the thermal averages [Eq. (26)] of colorless operators 
O the unitary matrix U{9) drops out. Since the density 
matrix Vg does not depend on the coset degrees of free- 
dom 9, the corresponding integral can then be trivially 
carried out 



S2 



dfi{e) ^ 4tt, 



(34) 



(27) and we obtain for the projected thermal averages 



{0)=- I dp{9) Z{6) {0)g , Z= I dp{e) Z{9), (35) 



where 



{0)e 



z{e) 



tr(PeO), Z{e) - TvVe (36) 



denotes the thermal expectation value for a fixed color 
angle 9. 

Furthermore, it is also convenient to use the basis in 
color space in which (Ta)"*" = e"^^^ is diagonal. In this 
basis we have 



Q=0,±1 '' 

and the density operator Vg [Eq. (32)] becomes 



Vg = cxp 
where 



-f3j2jdke-{k,e)af{k)ank) 



/3e"(k,0) ^ f3n{k)-i9a. 



(37) 



(38) 



(39) 



Since Vg is the (exponent of a) single particle operator, 
the thermal expectation values (. . .)g can be evaluated 
using Wick's theorem. In the standard fashion one finds 
for the partition function 



Z{0) = exp 



b'-^I 



dfcln[l + n„(k,0)] 



(40) 



where V is the volume of ordinary space and 2 — tu (k) is 
the number of independent polarization degrees of free- 
dom in three dimensions. Furthermore, 



n„(k,0) = (c'3^"(k,9) _ 1 



(41) 



are the finite-temperature Bose occupation numbers. 
The basic contraction is obtained as 

(af^(k)a^(q)) = (5"^i,,(k) (2^)35(k - q) n„(k, 9). (42) 



Expressing the gauge field in terms of the creation and 
annihilation operators one finds from Eq. (42) for the 
gluon propagator 



l + 2na(p,( 



{Afip)A^i<l)) = <5"^ t.,(p) i2nf6ip + q) 

(43) 

With these relations it is straightforward to calculate the 
thermal expectation value of the Hamiltonian using the 
same approximation as at zero temperature in Ref. [40], 
i.e., assuming a bare ghost-gluon vertex and calculating 
the energy up to two loops. 

To work out the effect of the color projection on the 
energy, let us for the moment ignore the Faddeev-Popov 
determinant in the Hamiltonian. We will later fully in- 
clude J [A]. Using the explicit form of the thermal gluon 
propagator, Eq. (43), and the same approximation as in 
Ref. [40] but putting J [A] = 1, one finds for the various 
pieces of the energy 



{HK)e = ^Jdqoj{q)[3 + 2Y,M^)], (44a) 

(^B)« = ||dg-^[3 + 25]n„(q)] 

V^/d,d/-(^-P)' 



16 J "■"" w(q)a;(p) 
3 + 2^(nQ(p) +na(q)) + 2 y^ Wg (p) n/j (q) 



a.fS 



^na(p)(n„(q)+n_a(q)) 



, (44b) 



{Hc)e 



g^N, 



V Jaqap[l + ici-pf]F{q-p) 

{5S[3 + 2^(n„(q)+n„(p))+2(^n„(q))(^n^(p) 

-^na(q)(n„(p) +n_„(p)) - 3 + ^ na(p)(n„(q) -n-a(q)) 

+ 13, V2- F(0)y dpdg^n„(p)(n„(q) - n_„(q)). (44c) 



I 

To simplify the notation, we have omitted the 0-depen- studies [37, 38] to have the infrared behavior 

dcncc of the occupation numbers Eq. (41). In Eq. (44c) , /,4 

-T (k — > 0) ^ 1/fc . (46) 

Accordingly, the integrand in Eq. (44c) becomes diver- 
is the non-Abelian color Coulomb potential. This quan- gent for p = q. Furthermore, the last term in Eq. (44c) 
tity is known from the lattice and also from continuum is manifestly divergent. However, one easily shows that 



these divergent terms disappear after color projection. 
For this purpose, we note that if one replaces in the 
Coulomb Hamiltonian He [Eq. (2)] the Coulomb kernel 
Fa hy 



-<(x,y)^<5"^<5(x~y), 



(47) 



the Coulomb Hamiltonian becomes the square of the total 
charge 



He -^ Q'Q" 



(48) 



This equivalence holds even when the Faddeev- Popov de- 
terminant is included, since J[A], being invariant un- 
der global color rotations, commutes with the total color 
charge operator Q". In momentum space the replace- 
ment Eq. (47) corresponds to 



9 



F{k) -^ (27r)3,5(k). 



(49) 



It follows that the singular p = q contributions to the 
double integral in {Hc)e [Eq. (44c)] are proportional to 
{Q"'Q°')e- However, this quantity has to vanish after color 
projection 



{Q^Q^) = ^j Aii{e)z{e){Q^Q'') 



0. 



(50) 



Therefore, the singular contributions that occur from the 
p = q part of the integrand vanish after color projection. 
We can explicitly eliminate these singularities by replac- 
ing the Coulomb potential F(k) by 

F(k) ^ F(k) - F{fd)V-\2TTf5{\i) 

= F{\.){l-V-\2nf5{\.)). (51) 

This replacement will, in particular, remove the last term 
of Eq. (44c). The kernel F and thus the color- projected 
Coulomb energy (-ffc)e is invariant with respect to a shift 
of the Coulomb kernel by a constant 

F(x,y)^^(x,y)+C. (52) 

This shift implies in momentum space 

F{\i) -^ F{k) + Ci2Trf6{k), (53) 

which obviously leaves F(k) [Eq. (51)] unchanged. 

B. Color projection in the thermodynamic limit 

The partition function Eq. (40) depends via the finite- 
temperature occupation numbers na(k) [Eq. (41)] on the 
color angle 9. The ^-dependence can be explicitly sepa- 
rated yielding 



z{0)^z{o)cM-v.ne)), 



(54) 



where 

Z{0)^cxpl2{N^ ^1)V f dk\n[l + n{k)]\ (55) 

is the partition function for vanishing "external" color 
field (9 = 0), with 

n(k) = n„(k)|e=o - n„^o(k) - (c^^^C^' - l) "' (56) 

being the corresponding thermal occupation numbers. 
The 0-dependence is entirely contained in the exponent 
of Eq. (54) , which is given by 

f{e) = 2 /"dfcln[l + 2(l-cos0)n(k)(l + n(k))]. (57) 

Note that the partition function is an even function in 
6. This property holds for the expectation value {0)e of 
any color singlet operator O. 

Consider now the total partition function Z, Eq. (35). 
In the integration domain 9 € [— 7r,7r] the function f{9) 
[Eq. (57)] takes its minimum at 9 ~ 0, where it vanishes 



f{9 = 0) = 0. 



(58) 



Due to the presence of the volume factor V, in the ther- 
modynamic limit T^ — > oo only small 9 values contribute 
to the integral Eq. (57). Therefore it suffices to expand 
the function f{9) to leading order in 9 yielding'^ 



Z{9) =Z{0)c- 



\ce' 



(59) 



where 



C = V,f"{Q), /"(O) =2 / dfcn(k)(l + n(k)). (60) 



With this representation for Z{9), the thermal expecta- 
tion value Eq. (35) of an observable O becomes 



(0) = ^ /dM(^)(0),e-^^«^ 



with 



Z = Z(0) / d^(e') c 



(61) 



(62) 



Due to the presence of the Gaussian, only small 9 val- 
ues contribute significantly to the integrals. Therefore it 
suffices to expand {0)q up to leading order in 9 



(O)e = (0)e=o + 0('^e' + --- 



Defining 



Z(0) 



dii{9) Z{9) 



Q2(n-1) 



(63) 



(64) 



^ The same expansion was used in Ref. [46] for the quark partition 
function. 



and using 



we obtain 



Z = Z(0)/i , 



(0) = (0)0=0 + 7^0^"^ 



(65) 



(66) 



Along the same lines, we can also expand the integration 
measure d/i(0) [Eq. (33)] to leading order in d and put 
the upper integration limit to oo. This yields for the 
integrals Eq. (64) 



1 

2^ 



d0 02ng-iCe^ 



1 (2n-l)!! 



V8^ C"" 



(67) 



Since l2li\ ^ 1/^ in the thermodynamic limit V — > oo, 
the second term in Eq. (66) can be omitted and we find 



(0)-(0)e=o. 



(68) 



This shows that in leading order in the thermodynamic 
limit the effect of the color projection can be ignored. In 
the following we will skip the subscript 6* = and (O) 
means (O)e=0; which is the unprojected thermal average. 
To include the Faddeev-Popov determinant we use the 
representation [47] 



J[A] = exp (-IJ^k A{-k)x{k)A{k) 



(69) 



where 



x(p) = ^ / d"P (i - (p • q) 



rf(p - q) rf(q) 

(p-q)' 



(70) 



is the ghost loop (curvature) and d{p) is the ghost form 
factor defined by 



((-^a)-)^i^(-^) 



(71) 



9 (-A) ' 
which satisfies the following Dyson-Schwinger equation 



d-\p) = --iAp), 



/,(p) = ^|d,[l-(p.q)2] 



2. d(p - q) 1 + 2n(q) 



(p - q)^ '^(q) 



(72) 



where a bare ghost-gluon vertex has been assumed. This 
equation differs from the zero-temperature case only by 
the replacement of the gluon propagator by its finite 
temperature counterpart Eq. (43). The representation 
Eq. (69) is valid up to two loops in the energy, which is 
the order considered in the present paper. 

With the inclusion of the Faddeev-Popov determinant, 
the thermal expectation value of the Hamiltonian (Hym) 
given by Eqs. (44) simplifies ioi 9 = to 



{Hym) 

where 


= {n! - 


1) 


■2-V e, 


e 


= eK + eB + ec, (73) 


ex = 


i/.,. 


([" 


^^) + r 


^(q)] 


1 + 2n(q) 
w(q) 


-2x(q)} 

(74a) 



1 f 2 i + Mq) , g^N, 

es = - / dg q — — + — — 

4 ./ oj(q) 32 



J dp dq[3 - {p ■ qf] 



1 + 2n(p) 1 + 2n(q) 



5^ 
16 



(x)(q) 32 7 "" ^'- \^ ^/ J a;(p) u;(q) 

I dp dq[l + (p . q)2] ^(^^{ [^'(P) + X'(P) - X(p)x(q)] (1 + 2n(p)) (l + 2n(q)) 



c^(p) c^(q) + 2x(p)[c^(q)(l + 2n(p)) - c^(p)((l + 2n(q))] } 



(74b) 



(74c) 



are the energy densities per degree of freedom. 



IV. FINITE TEMPERATURE VARIATIONAL 
PRINCIPLE 

Our ansatz for the density operator [Eqs. (20) and (27)] 
contains a so far arbitrary kernel r2(k), which we deter- 
mine now by minimizing the free energy J^, Eq. (21). 
Instead of varying T with respect to r2(k), it is more 
convenient to take the variation with respect to the finite- 



I 

temperature occupation number n(k) [Eq. (56)], which is 
a monotonic function of 51 (k) for 51 (k) > 0. Variation of 
T with respect to n(k) yields 



^(k) = 



(5e[n] 
Sn{k) 



(75) 



which identifies r2(k) as the quasi-gluon energy. 

So far, the kernel w(k), which defines the vacuum wave 
functional Eq. (8) and thus our basis of the Fock space, is 
completely arbitrary and, in principle, we could use any 
positive-definite kernel w(k). As long as we include the 



complete set of states and do not introduce any approxi- 
mation, the thermal expectation values will be indepen- 
dent of u;(k). However, due to approximations necessary 
as, for example, the restriction to two loops, the thermal 
averages will depend on the u;(k) chosen and the optimal 
choice is obtained by extremizing the free energy Eq. (21) 
with respect to cj(k) 



5F 



= 0, 



<5cj(k) 
which yields the finite temperature gap equation 



;(k) 



where 

/(o) 

/(k) 



g^N, 



g'N, 



x2(k) + /W+/(k), 



<I9 ^ ^^''„;^^' [1 + 2n(q)] 



(76) 



(77) 



u;(q) 
dgF(k-q) 



{ 



1 + (k ■ q)^ 

u;(q) 



[1 + 2n(q)] [c.2(q) - ^.^(k) + (x(q) - x(k)) 
-2c.(q)[x(q)-x(k)]}. 



(78) 



These loop integrals, as well as the ones of the curvature 
[Eq. (70)] and of the ghost Dyson-Schwinger equation 
[Eq. (72)], are ultraviolet divergent and need to be regu- 
larized and eventually renormalized. 

Inserting the explicit expressions Eqs. (74) for the en- 
ergy densities e into Eq. (75) and using the gap equation 
(71), one finds 






Uk) = 



SqFik-q] 



1 + (ic ■ q)^ 
uj{q) 



To carry out the renormalization we will have to deal 
with both finite- and zero-temperature solutions. To 
avoid confusion, in the following we will explicitly in- 
dicate the temperature dependence by writing a;(k, T), 
d(k, T), ... instead of w(k), d(k), . . . 

At very large momenta |k| ^ T the temperature 
should become irrelevant. Indeed, the temperature de- 
pendence of the loop integrals (which is due to the finite- 
temperature occupation numbers n(k)) does not give 
rise to additional UV singularities. Therefore the zero- 
temperature counterterms are sufficient to eliminate all 
UV singularities. 

Adding the zero-temperature counterterms, see 
Ref. [41], and carrying out the renormalization as de- 
scribed in Ref. [48], one arrives at the following renor- 
malized gap equation 



^'(k,T) 



k^ + x' (k, T) + A/(2) (k, T) + Co + /(°) (T) 
-f2x(k,T)[A/«(k,r) + ci]+J(k,r), 

where we have introduced the abbreviations 



x(k,r) 

/(0)(r) 

A/(')(k,r) 



x(k,r) 

/(0)(r). 



/(°'(r = o), 



(81) 



/W(k,T)-/«(A.^,r = o), 



[l + 2ri(q)] . and defined the following loop integrals 



i^ 



1 + (k • q 



/ W (k, T) = ^ y dg F(k - q, T) i±^^ { Mq, T) - x(q, T)] ' - Hk, T) - x(k, T)] ' } , 



/(k,T)^^/d,F(k-q,T)^ + (^-^)^ 



^(q,T) 



2n(q) {c.2(q, T) - c.2(k, T) -|- [x(q, T) - x(k, T)] '} 



(82) 



(83) 



In addition, the renormalized equation (72) for the 
ghost form factor reads 



1 



d(k,T) 



1 



di^id,T = 0) 



+ /<i(/Xd,r = 0) -/d(k,T). 

(84) 

The renormalized Dyson-Schwinger equations (80) and 
(84) contain the finite renormalization scales, fii ~ 
Mx' /^"' t^d, and the renormalization constants g{p-i), 
x(Mi)) co(/Xi), and Ci(/Xi). The last two originate from 
the counterterms in the Hamiltonian and xip-i) from the 



r 



renormalization of the Faddeev-Popov determinant. In 
particular, in Ref. [48] it was shown that for ij-uj — fJ-d — 
the value Ci = is required in order that the 't Hooft loop 
obeys a perimeter law and is also favored by the varia- 
tional principle. It was also found that the parameter 
Co has no influence on the IR- or UV-behavior of the re- 
sulting solutions and influences only the mid-momentum 
regime of w(k). The choice of renormalization conditions 
for our study at finite temperature will be discussed in 
Sec. V. 



The Gribov-Zwanziger confinement scenario assumes 
that d~^{0,T = 0) = 0. For practical reasons, in the 
present paper we will assume a small but finite d~^{0, T), 
which results in a massive gluon propagator, referred to 
as subcritical solution in Ref. [41]. This solution does 
not provide a confinig Coulomb potential, but for phe- 
nomenological purposes may be as useful as the critical 
confining solution d^^{0) = (see Ref. [41] for further 
discussions). One can give arguments that a (i~^(0) ^ 
is the result of an improper treatment of the Gribov prob- 
lem [33, 49]. In fact, it was explicitly demonstrated in 
1 + 1 dimensions [50], and also arguments were given 
for 3 + 1 lattice gauge theory in Landau gauge [49], 
that extending the functional integral over the trans- 
verse gauge field to higher Gribov regions reduces the 
infrared strength of the ghost form factor, pushing d~'^{0) 
to higher values. Based on this observation, it was ar- 
gued in Refs. [33] and [49] that choosing different values 
of (i^^(O) corresponds to different "gauge fixings" . (After 
all, a complete gauge fixing implies also the restriction to 
the fundamental modular region, which is a subset of the 
first Gribov region.) Presumably, in more than 1 + 1 
dimensions the restriction to the fundamental modular 
region requires d~^{0) — 0. In any case this value is 
required for a linearly rising Coulomb potential, which 
is a necessary condition for confinement in the Gribov- 
Zwanziger confinement scenario [42, 43]. Thus, if d^^(O) 
is kept finite for technical reasons, it has to be kept small 
to stay close to the physical confining limit d'^{0) — 0. 



V. NUMERICAL RESULTS 

As shown in Sec. Ill A, the color projection removes 
the zero mode from the Coulomb potential, see Eq. (51). 
In the continuum it is replaced by 



100 



g^F{k) = lim 



d2(k) 



«^^0 k^ + £2 



[l-exp(-kVe')], 



(85) 



were we used the approximation [cf. Eq. (6)] [40] 



«(x,y)) ^^ /d3z((x,al(-DV)-ilc,z)) 

xi-V\{{z,c\{--DV)-'\y,b)), (86) 

with the external (• • ■ ) referring to the thermal average. 
In the confining limit 



g'Fik) = 



d\k) 



k^O 



87rcrc 



(87) 



and the single, quasi-gluon energy in Eq. (79) is infinite 
at all temperatures, which is certainly an artifact of our 
approximation, since at least for large ]k] the quasi-gluon 
energy should be finite due to asymptotic freedom. For 
infinite J7(k) the finite-temperature occupation numbers 
n(k) [Eq. (56)] vanish at all temperatures and there is no 



3 




FIG. 1. Low momentum (IR) behavior for j3 = 0.5, e — 0.01, 
and dg^ = 0.0005, 0.0040 of solutions for gap equation for cu. 
In the limit e — > and/or /? ^ c« the solutions do not change 
qualitatively. The critical solution corresponds to d^^ = 0, 
and the solution with d. 







0.0005 is close to critical. The IR 
limit is weakened as dg^ increases and the ghost propagator 
becomes massive. 



T3 




0.01 - 



FIG. 2. Same as Fig. 1 for the ghost form factor d. 



finite-temperature phase transition. Thus the presently 
used approximations are inappropriate for the strictly 
confining solution. As discussed in Ref. [41], without 
the approximation of Eq. (86) there are no strictly con- 
fining solutions in the sense of Eq. (87) in the variational 
approximation, and in this case the r2(k) are finite (for 
finite k) and at finite temperature a non-trivial solution 
with n(k) ^ is expected. In this case we solve the set 
of finite-temperature Dyson- Schwinger equations numer- 
ically on a momentum grid. The non-confining solutions 
depend on the renormalized coupling or, alternatively, 
the value of d(k = 0,T) [cf. Eq. (72)]. As discussed in 
Sec. IV, the Dyson-Schwinger equations are renormal- 
ized by subtraction at zero temperature to account for 
temperature-independent counterterms. In the numer- 
ical computation, however, it is very difficult to solve 
these equations at fixed T unless subtracted at the same 
value of T. Thus, in the numerical results that follow all 
subtractions will be done at finite T. In particular, when 



solving for (i(k, T) [see Eq. (72)] we use 



d(k,T) 



d{^icuT) 



+ Id{lJ.d,T) 



IdiKT). 



Comparing this with Eq. (84) we have 



1 



1 



d{fid,T) d{fid,0) 



+ Id{pd,0)-Id(pd,T), (89) 



with fj,d chosen to be the lowest point on the momentum 
grid, which corresponds to fj,d = in the infinite volume 
limit. In other words, we fix dg with the temperature: 
At each temperature we thus control the distance to the 
confining limit of the color Coulomb potential. This im- 
plies that the mass scale which enters IdifJ'd ^ 0,T) on 
the r.h.s. of Eq. (89) depends on T. Similarly, the numeri- 
cal stability of the solution of the gap equation for w(k, T) 
[Eq. (80)] requires that we use temperature-dependent 
renormalization constants, i.e., in Eq. (81) instead of 
subtracting at T = we subtract at finite T, so that 
/i^ = fJ-xiT) and /z^j = /Z(^(T). In particular, we use a sin- 
gle renormalization scale and set ^(T) = fJ.io{T) = Mx(^)- 
This implies that we renormalize the gap equation at a 
finite momentum fi^j 7^ 0. By renormalizing at fi^j — 0, 
one would be enforcing a particular IR limit of the so- 
lution of the gap equation, which could turn out to be 
incompatible with the finite-temperature equation. In- 
stead, by choosing /i^j away from the IR limit the value 
obtained from solving numerically for w(0) = w(k = 0, T) 
will serve to illustrate the onset of a phase transition. 

To search for the phase transition we proceed as fol- 
lows. We start with a small but finite IR regulator e [see 
Eq. (85)], and at given, small T (large /3) we choose a 
solution close to a critical one. A typical case is shown 
in Figs. 1 and 2. In all figures, physical, dimensional 
quantities are plotted in units of /x(T). With fixed e and 
do we increase T (decrease /3) and study both lu and d. 
A series of computations of w(0) as function of /? for 
e = 0.01 and d^^ in the range [0.0005-0.0060] is shown 
in Fig. 3. As dg^ increases, the solution becomes less 
critical, i.e., less IR enhanced, and the Coulomb poten- 
tial moves away from the confining limit of Eq. (87). In 
this phase, as shown in Fig. 3, there is an abrupt change 
in the gap function w(0) at a critical temperature which 
decreases as the solution become weaker in the IR. Please 
note that the results shown for different values of d~^{0) 
correspond to different physical scales. We also studied 
the dependence on e, as the limit e ^- should be taken 
to approach the infinite volume. Starting from a massive 
solution at zero temperature, e.g., with dQ = 0.060, we 
increase the temperature and decrease e. In particular, 
we solve the equations for e = 0.0001 and e — 0. The lat- 
ter choice is possible, since with an IR finite ghost dress- 
ing function (d^^ ^ 0) the Coulomb potential in Eq. (87) 
does not lead to an IR singularity in the integrals. The 
results are shown in Figs. 4 and 5. Even though, in nu- 
merical simulations with a finite momentum grid, one 
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FIG. 3. ijj{0) as a function of temperature for e = 0.01 and 
dg^ = {0.0060, 0.0040, 0.0029, 0.0015, 0.0005}. The phase 
transition is clearly visible and becomes stronger and moves to 
lower temperatures as d^^ increases. For d^^ ^- the phase 
transition disappears, i.e., moves to infinite temperatures. 
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FIG. 4. cj(0) as a function of temperature for d^ = 0.0060 
massive solution for e = 0.01 (as in Fig. 3) compared with 
with solutions for e — 0.0001 and e = 0. 
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FIG. 5. Same as in Fig. 4 for d^^ = 0.0029. 



cannot reach the exact critical limit of e = 0, dn 
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FIG. 6. Gap function cj(k) as a function of temperature 
(dp ^ = 0.0029, e = 0). As the temperature increases u}(k) 
becomes less IR enhanced. 



(which implies i^(k) ex 1/k^ at zero temperature), the 
numerical results shown in these figures are consistent 
with the anticipated (see Sec. I) disappearance of the 
phase transition in the Coulomb phase. In particular, 
we observe that, as d^^ decreases, the gap function w(0) 
grows and eventually becomes infinite as it is the case 
in the zero-temperature limit. In other words, the phase 
transition moves to infinite temperature (/? — > 0) as the 
gluon self-energy becomes infinite. In Fig. 6 we show the 
evolution with temperature of a;(k) as a function of gluon 
momentum. As expected, the gap function becomes less 
IR enhanced as temperature increases {(3 decreases). 



Strictly speaking, the variational Coulomb gauge 
model which leads to d^ ^ is not confining and thus 
only loosely related to QCD. The Gribov-Zwanzinger 



confinement scenario is reached in the dr, 



(and 



e — > 0) limit. We find it amusing, however, that the 
quasi-gluons which for d^ y^ are deconfined at all tem- 
peratures behave similar to the physical gluons both be- 
low and above Tc- We have found that quasi-particle, 
gluonic excitations built on top of such a subcritical vac- 
uum lead to a sharp transition in the above correlations 
functions. To solve the Dyson-Schwinger equations we 
used temperature dependent renormalization conditions. 
This results in the phase transition point in Fig. 3 moving 
with dQ^. By fixing the critical temperature to a physical 
value, if known, this variation could be used to determine 
the function fi{T) and, ultimately, the temperature de- 
pendence of these correlation functions. As one tunes the 
zero-temperature solution to approach the critical limit, 
dp — >■ 0, the phase transition moves to infinite temper- 
atures. Even though thermal excitations are restricted 
to color single states, contribution to the partition func- 
tion from two gluons is 0{1/V) compared to that of a 
glueball. The former are thus expected to make negli- 
gible contribution in the thermodynamical limit, which 
explains why a confining Coulomb potential at zero tem- 
perature remains confining at finite temperatures [37]. 
However, the confining potential can bind gluons into 
color-singlet glueballs and a phase transition could be 
observed, for example in a change of the radius of the 
glueball wave function. We will consider such a mixed 
glueball/quasi-gluon phase in the forthcoming work. 



VI. SUMMARY AND OUTLOOK 

We studied the temperature dependence of QCD cor- 
relation functions with a variational ansatz for the gluon 
density matrix in the Coulomb gauge. The resulting one- 
loop Dyson-Schwinger equations for the gluon propagator 
and ghost form factor (i(k) were solved numerically, as- 
suming a subcritical behavior, i.e., d^^{0) ^ 0, which, 
however, was chosen close to the critical one, d^^{0) = 0. 
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